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ABSTRACT 

We address the fundamental question of matching the rest-frame if -band lumi- 
nosity function (LF) of galaxies over the Hubble time using semi-analytic models, after 
modification of the stellar population modelling. We include the Maraston evolution- 
ary synthesis models, that feature a higher contribution by the Thermally Pulsating - 
Asymptotic Giant Branch (TP-AGB) stellar phase, into three different semi-analytic 
models, namely the De Lucia and Blaizot version of the Munich model, MORGANA and 
the Menci model. We leave all other input physics and parameters unchanged. 

We find that the modification of the stellar population emission can solve the 
mismatch between models and the observed rest-frame A-band luminosity from the 
brightest galaxies derived from UKIDSS data at high redshift. For all explored semi- 
analytic models this holds at the redshifts - between 2 and 3 - where the discrepancy 
was recently pointed out. The reason for the success is that at these cosmic epochs 
the model galaxies have the right age (~1 Gyr) to contain a well-developed TP-AGB 
phase which makes them redder without the need of changing their mass or age. 

At the same time, the known overestimation of the faint end is enhanced in the 
A-band when including the TP-AGB contribution. At lower redshifts (z < 2) some of 
the explored models deviate from the data. This is due to too short merging timescales 
and inefficient 'radio-mode' AGN feedback. Our results show that a strong evolution in 
mass predicted by hierarchical models is compatible with no evolution on the bright- 
end of the A-band LF from z=3 to the local universe. This means that, at high redshifts 
and contrary to what is commonly accepted, A-band emission is not necessarily a good 
tracer of galaxy mass. 

Key words: methods: numerical - methods: statistical - galaxies: formation - galax- 
ies: evolution - stars: AGB 



1 INTRODUCTION 

Models for the formation and evolution of galaxies in a 
Cold Dark Matter universe (e.g. the so-called semi-analytic 
models) predict the intrinsic properties of galaxies, such 
as ages, metallicities, stellar masses, star formation rates, 
etc., after having tuned a number of free parameters 
that make u p for the p oorly known aspects of baryonic 
physics (see iBaughl 120061 . for an extensive review). The 
comparison between models and observations helps con- 
straining these parameters and robust statistical tools have 
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been r ecently used to a c hieve t his goal dKampakoglou et al.l 



| 2008l: IHenriques et al.1 120091: iHenriques fc Thomas! l2oiof 



bower et al.ll20ld : lLu et al.ll2010h 



The results of these comparisons are very sensitive to 
the spectro- modelling of the stellar component. Either for 
extracting galaxy properties such as mass, age, star for- 
mation rate from data, and compare them to the intrin- 
sic quantities of the semi-analytic models, or for calculating 
the spectra of semi-analytic galaxies and compare it with 
the observed light, the details on how the stellar modelling 
is performed influence the final result. 

In order to obtain the spectral energy distribu- 
tion (SED), or specific broad-band luminosities, of a 



© 2009 RAS 



2 Henriques et al. 



model galaxy of given mass and star formation his- 
tory, evolutio n ary p o pulations s y nthes is (EPS) models 
(e.g. iTinslevI Il97l Bruzual A.I Il983t lBu"zzonil Il989t 
Bruzual A. fc Charlotl 1 19931. IWorthevlll994 iVazdekis et al 



1996 



Fioc fc Rocca-VolmerangeT 



1997 



Maraston 



19981. iLeitherer et all 1 19991 IBruzual fc Charlotl I2003L 



Thomas et al.1 120031 . IMarastonl 120051 (M05), IConrov et all 



20091 ) are adopted. By relying on stellar evolution theory 



and model atmosphere calculations or empirical libraries, 
EPS models provide the expected spectral energy dis- 
tribution of a galaxy of given mass and star formation 
history. 

In galaxy formation models the unit single burst EPS 
models or Simple Stellar Populations (SSPs) are used to 
model coeval stars with a homogeneous metallicity, after 
adopting an Initial Mass Function (IMF). The total stellar 
emission from the synthetic galaxy composite population is 
then obtained by combining SSPs. Hence, what matters on 
the final result in terms of stellar evolution are the properties 
of the simple stellar population models. 

It is clear then that this modelling is a crucial aspect 
of galaxy formation and evolution theory. Uncertainties in 
the conversion between masses/ages and light can create 
artificial discrepancies, which in turn could drive into dif- 
ficult attempts to modify the parametrization of the com- 
plicated physics of gas cooling, star formation or feedback 
to account for this mismatch. The approach we take i n this 
paper, following our previo us work (|Tonini et al . 2009, 2010; 
iFontanot fc Monacdl201Cl ). is to check the impact of modi- 
fying the input stellar population model in galaxy formation 
models following recent progress in the literature. 

At present, the highest source of discrepancy between 
different SSP models is the treatment of the Thermally- 
Pulsating Asym ptotic Giant Branch (TP - AGB) phase 
llMarastonl 1 19981 , M05, iMarigo et all 120081 . IConrov et ail 
12009] ). The TP-AGB phase is the last luminous phase in the 
Hertzsprung-Russel diagram before intermediate-mass stars 
evolve to their final destiny as planetary nebulae and white 
dwarfs. TP-AGB stars are very luminous and cool. Their 
emission affects the integrated model spectra at wavelengths 
larger than ~ 6000A, peaking around the J, H, K bands 
(M05) and a recent study has hig hlighted their importanc e 
also long-ward the rest-frame K (iKelson fc Holder] [2010). 
Due to difficulties in the stellar modelling of this phase, 
whic h in turn are due to m ass-loss and the pulsating regime 
(see llben fc Renzinil 1 19831 . for a review) , generally stellar 
tracks did not include the full TP-AGB, so did not stel- 
lar populat ion models based on these tracks (see M05 for 
discussion). IMarastonl l| 19981 ) and M05 include the TP-AGB 
semi-empirically calibrating the theoretical energetics with 
Magellanic Cloud clusters, an approach n ow adopted for in- 
cluding the TP-AGB phase in isochrones (IMarigo fc Girardil 
l2007l ; lBruzuaill2007l) . 

The galaxy formation models we use in this study 
do not include, in their standard stellar populations, the 
TP- AGB stellar phase. As pointed out by several pa- 
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) the inclusion of TP-AGB stars provides 



of this inclusion on model predictions, one needs to consider 
a statistically significant sample of model galaxies, covering 
both a wide range of K-band luminosities and redshifts. This 
test is particularly important when the photometric proper- 
ties of high-z (i.e. 2<z<3) galaxies are considered, since we 
expect them to be dominated by young stellar populations. 

Semi-analytical techniques represent an obvious tool 
to perform such test. Predictions from these models have 
already been compared with the evolutio n of the ob- 
serve d jf-band Lumi n osity F unction (LF) (iPozzetti et al.l 



120031; ICimatti et al.1 |2004| ; iKitzbichler fc White! |2~ 



Cirasuolo et al.ll2010l ). These works consistently found a lack 
of bright sources at high redshift. This apparent mismatch 
is being referred to as one of the strongest discrepancies be- 
tween models and data (in particular in con nection with the 
evolu tion of the stellar mass function, see IFontanot et al.l 
2009b for a critical review of the latter issue). However, 
these comparisons involved spectro-photometric codes based 
on stellar tracks where the effect of TP-AGB stars was not 

taken into account. 

T he firs t test of this kind is performed in lTonini et al.l 

20091. |2010|). They run the GALICS semi-analytic model 



Hatton et all [2003) using M05 as input EPS and show 



an enhancement of the near-infrared emission of galaxies 
dominated by ~1 Gyr old populations. To test the influence 



that the optical-to-near-IR colours of z~2 galaxies can be 
matched by this type of models, with the original GAL- 
ICS, which was based on a pre-TP-AGB EPS, failing 
to match the observation s by a large margin. Moreover, 
IFontanot fc Monaco! (|2010l ) introduced the M05 stellar pop- 
ulations in MORGANA obtaining a good match on the number 
density of extremely red objects (EROs) at high redshift. 

Here we investigate whether the inclusion of the TP- 
AGB phase has also an impact on the inability of semi- 
analytic models to match the galaxy rest-frame if -band 
luminosity function at high-reds hift, namely the UKIDSS 
data from IcTrasuolo et al.l l|2010l) . This homogeneous data 
set covers an ideal redshift range (0<z<3) for this test. 

To perform this analysis, we use three differ- 
ent semi-analytic mode ls of galaxy formation: the 
|Pe Lucia fc BlaizotJ (120071) version of the Munich model, 
morgana l|Monaco et al.ll2007l ) and the lMenci et al.l (|2006l ) 
semi-analytic model. We compare the predictions obtained 
for the properties of the galaxy population using both stel- 
lar populations with and without a full treatment of the TP- 
AGB phase. For the three models (|Pe Lucia fc Blaizotll2007l ; 
iMonaco et al.ll2007l ; lMenci et a l. 2006) ou tputs a re produced 
using, resp ectively IBruzual fc Charlotl l|2003l ). ISilva et ail 
i|l998l) and IBruzual A. fc Charlotl (|l993l ). and M05. 

This paper is organized as follows. In Section [2] we 
briefly describe the semi-analytic models used in this study, 
and we explain how the M05 EPS is implemented in each 
galaxy formation model. In Section [3] we describe the data 
used for comparison and clarify where the impact of the M05 
models is expected to be found. Section|3]presents the results 
for the evolution of the rest-frame near-infrared luminosity 
function and in Section [5] we summarize our conclusions. 



2 THE SEMI-ANALYTIC MODELS 

A clear advantage of the work presented in this paper is 
the use of three semi-analytic models developed by inde- 
pendent groups and implementing different techniques for 
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the description of the physics controlling galaxy formation 
and evolution. This allows us not only to assess the im- 
pact of the TP-AGB phase, but also to understand the in- 
terplay between the new ingredient and the other assump- 
tions regarding galaxy physics on the light of the same 
stellar evolution ba c kgroun d. In particular we consider the 
|Pe Lucia fc BlaizotJ (|2007f ) version of the Munich model; 
MORGANA (o riginally described in iMonaco et~aH 12007 
and u pdated bv lLo Faro et al1l2009l ): and the iMenci et al. 
(|2006h model. 

The backbone for all models is a description of the red- 
shift evolution of the mass and number density of dark mat- 
ter halos in terms of their merger history (the so-called 
merger trees). The evolution of the baryonic component 
hosted by these halos is then followed by means of an ap- 
proximated set of simplified formulae, aimed at describing 
the physical processes acting on the gas (such as gas cooling, 
star formation and feedback) in terms of the physical prop- 
erties of each model galaxy and/or its components (i.e the 
stellar, hot and cold gas content and distribution). These 
analytical recipes include a set of parameters which are usu- 
ally calibrated against a well defined subset of low-redshift 
observations. 

The three models adopt different techniques to describe 
the dark matter merger treef] and slightly different cosmolo- 
gies. However, we do not expect these to have a sig nificant 
effect on our conclusions (see e.g. IWang et al 1 I2008T ). 

On the other hand, the different star formation histo- 
ries and the corresponding distribution of ages and the mass 
build up in the models do matter. In the following, we will 
briefly account for differences between the models, focusing 
in particular on the AGN feedback and the merging time 
scales, the processes most relevant for the evolution of the 
bright-end of the jf-band LF. For more details on the treat- 
ment of these physical processes in the diffe rent models we 
refer t he reader to the original papers, and to lDe Lucia et aD 
l|2010l ); Fontanot et al. (2009b) for recent comparisons. 



2.1 AGN Feedback 

The recipe adopted to describe AGN feedback is of crucial 
importance, since it largely determines the stellar population 
properties of the most massive galaxies, whos e evolution is 
the f ocus of our paper. Recent studies (e.g. ICroton et al.l 
2006) assume that the growth of Super-Massive Black-Holes 
(SMBHs) at the center of model galaxies follows two chan- 
nels, a "bright-mode" (or "quasar-mode") and a "radio- 
mode", related to the efficient production of radio jets. The 
"quasar-mode" is fueled by merger driven instabilities, it is 
the dominant channel in terms of black hole growth and 
can be effective in producing feedback at early times (where 
merging rates are high). The "radio- mode" is less impor- 
tant in terms of SMBH growth but is responsible for star 
formation quenching at low redshift. 



1 The Munich model uses merger trees extracted from a di- 
rect N-Body simu lation of a cosmolog ical volume (the Millen- 
nium Simulation, jSpringel et alj|2005l) ). MORGANA uses the La- 
grangian semi-ana lytic code PINOCCHIO {Monaco et al.ll2002h and 
IMenci et afl (120061 1 uses Monte Carlo realizations of merger trees 
based on the halo merging probability given by the Extended 
Press-Schechter formalism. 



The details of the implementations of the two modes 
differ between the m odels considered in this study (see e.g. 
iFontanot et al.ll2010b for a de t ailed discussion about mor- 
gana and lDe Lucia fc BlaizotJ (|2007f )y While the net effect 
of the "quasar-mode" is quite similar between the various 
models, even if the implementations are slightly different, it 
is the "radio- mode" , to be mostly responsible for differences 
in the galaxy stellar populations towards low redshift. 

In the Munich model, the "radio-mode" feedback is 
the result of quiesc ent gas accretion from a static hot halo 
jCroton et al.l2006l ). with no triggering mechanism required. 
In MORGANA, the "radio-mode" is due to the accretion (at 
very low rates) of col d gas from a reservoi r surrounding 
the central SMBH (see IFontanot et all (|2006t ) for more de- 
tails). Note that some amount of star formation is required 
to destabilize the gas in the reservoir. Hence, star forma- 
tion is not completely quenched. This residual star forma- 
tion causes galaxies to have colours that are too blue at low 
redshift with respe ct to both observations and other models 
|Kimm et al.ll2009h and contributes to an exces sive build up 
of ma ssive objects at later times. Finally, the IMenci et alj 
(2006) model does not include "radio- mode" feedback. For 
this reason, at low redshift, massive objects always have on- 
going star formation, which causes an excessive mass build 
up in these objects. Relevant to our work is that this results 
in an over-prediction of the bright tail of the K-band LF, as 
we will show in Section 2] 



2.2 Merging Times 

Dark matter substructures and their clustering have relevant 
consequences on the evolution of galaxies. Gravitational pro- 
cesses such as dynamical friction and tidal stripping affect 
the morphology, the stellar and the gaseous content of galax- 
ies. Two-body mergers are even more extreme processes, 
leading to the formation of a new object, whose final prop- 
erties depend on the properties of the progenitors. 

In the Munich model, dark matter substructures in the 
N-body simulation are explicitly tracked down until tidal 
truncation and stripping re duce their mass below the resolu- 
tion l imit of the simula ti on foe Lucia et al] [2004; Ga o et al] 
12004 ). In IMenci et all (120061 ) dark matter is followed us- 
ing a Press-Schechter formalism and satellite halos are par- 
tially disrupted as the density in their outer parts becomes 
lower than the d ensity of the host h alo within the pericentre 
of its orbit (see IMenci et al.1 (|2002t ) for details). After this 
point the merging time of the satellite in both m odels is 
computed using the classical IChand rasckharl (1 19451 ) dynam- 
ical friction approx i matio n. It is worth stressing that the 
|Pe Lucia fc BlaizotJ (|2007f ) model includes an additional pa- 
rameter in this formula, effectively doubling the expected 
merging times. This value was introduced to reduce the 
slight excess of bright galaxies that would be produced oth- 
erwise. MORGANA does not track explicitly dark matter sub- 
structures and assumes that satellite galaxies merge onto 
central galaxies after a dynamical friction time- scale which is 
compu ted using analytic formulae proposed bv lTaffoni et al.l 
(120031). 

|Pe Lucia et al.1 (|2010h compare different approxima- 
tions for the dynamical friction mergin g time-scales used in 
semi- a nalytics. They find that while the De Lucia fc BlaizotJ 
|2007l ) recipe is in good agreement with some recent re- 
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suits based on iV-body-simulations (|Bovl an-Kolchin et al.l 
2008),the lTaffoni et al] (|2003l ) formulae predict significantly 
shorter mer g ing t imes. Note that the same is true for 
iMenci et alj ll2006|) with merging times two times shorter 
than in |Pe Lucia fc Blaizotl (|2007r ) . 

Despite the overall agreement between different mod- 
els in terms of the mass build up found by Fontanot et 
al. (2009b), it can be seen that morgana shows an ex- 
cessi ve build up of ma ssive galaxies at late times. We ex- 
pect [MencT^£^l] (|2006l ) to show a similar behaviour. This is 
due to the combined effect from the enhanced merger activ- 
ity and the ongoing star formation due to inefficient AGN 
feedback at low redshift. This will affect our results with 
both models over-estimating the number density of bright 
A'-band objects at later times (see Section 2}. 

2.3 Implementation of the M05 models 

As recalled in the Introduction, the spectra of galaxies in 
semi-analytic models are obtained by means of spectro- 
photometric population synthesis models. The implemen- 
tation of the M05 models is straightforward in these semi- 
analytics. We use SSPs corresponding to four metallicities, 
1/20 Zq, 1/2 Zq, Zq and 2 Zq, which despite not being 
exactly the same as for the stellar populations previously 
used (since the input stellar tracks of the M05 models are 
different, see M05 for details), cover a similar range and 
are as coarse. Therefore, this difference has no impact on 
our predictions. The same IMF that was previously adopted 
in the va r ious s emi-analytic models is retained, namely the 
IChabrie j (|2003l) for MORGAN^, an d the Munich model, and 



Salpeterl (|l955h for lMenci et ail J2006). 



The predicted luminosities are then corrected for dust 
extinction. For all models we keep these prescriptions un- 
changed. The different treatment of dust extinction has non- 
negligible effects o n the predicted magn itudes and colors, es- 
pecially at z > 2 l|Fontanot et alj|2009h . However, since the 
rest-frame isT-band emission is relatively insensitive to dust 
attenuation, we do not expect these differences to substan- 
tially affect our results. 

Finally, It is worth stressing that we keep all other as- 
sumptions and parameters of the semi-analytic models as in 
their original formulation. Therefore, we can highlight any 
modification due just to the change in the stellar population 
libraries. 



3 A CHALLENGING ISSUE: THE OBSERVED 
REST-FRAME if-BAND LUMINOSITY 
FUNCTION AT HIGH REDSHIFT 

In this paper we focus on a well documented discrep- 
ancy between semi-analytic models and observations, the 
inability of the models in matching the observed red- 
shift evolution of the rest-frame near-infrared galaxy lu- 
minosity function dPozzetti et al.ll2003l ; ICimatti et al.ll2~004l ; 
iKitzbichler fc White! 120071 ) . This h as recently been con - 
firmed over a wide redshift range i|Cirasuolo et al.ll2010l ). 
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Figure 1. A simple stellar population spectral energy distribu- 
tion from M05 is compared with the equivale nt predictions from 
model s used in semi-analytic models (here IBruzual fc CharlotJ 
as an illustrative example). The plot refers to a 1 Gyr old 
population with solar metallicity. The TP-AGB stellar phase in 
the M05 models is significant between 0.2 and 1 Gyr. For similar 
plots and discussion see M05. 



This paper uses a data set from the Ultra Deep Survey 
(UDS), the deepest survey from the UKIRT Infra- Red Deep 
Sky Survey (UKIDSS), containing imaging in the J- and 
if-bands, with deep multi-wavelength coverage in BVRi'z' 
filters in most of the field. The sample contains ~ 50,000 
objects over an area of 0.7 square degrees, with high com- 
pleteness down to K ^ 23. ICirasuolo et al.l (2010) find that 
the space density of the most massive galaxies at high red- 
shifts (above 2) is under-predicted by semi-analytic models, 
in other words the theoretical luminosity function lacks the 

brig htest sources i n the near-I R. 

iTonini et all l|2009l , l20ld ); iFontanot fc Monacol (|2010l ) 
showed that the number of bright K-haxiA objects at high 
redshift in semi-analytic models can be increased by includ- 
ing the M05 models with their treatment of the TP-AGB 
phase of stellar evolution. We briefly recall here the origin of 
such an effect. The M05 models, predict that young popula- 
tions have a significant contribution to the near-infrared. For 
the luminosity function analysis, the differences are expected 
to be more significant at high redshift where a larger fraction 
of the galaxy population contains young stars. In Fig. [T] we 
plot the spectral energy distribution (SED) for a population 
with 1 Gyr of age and solar metallicity, using M05 and the 
stellar populations previously implemented in semi-analytic 
models (see M05 for similar plots and discussion). In the K- 
band the M05 model predicts more than twice the emission, 
hence affects the prediction of semi-analytic galaxies at high 
redshift. 

In concluding this section, some words of caution must 
be giv en on the data/model comparison. In lCirasuolo et all 
( 2010 ) the galaxies have photometric redshifts, which were 
obta ined by fitting empi r ical a s well as synthetic templates 
from IBruzual fc Charlotl (2003). For consistency, photomet- 
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ric redshifts and rest frame magnitudes should have been 
derived for the data using the same stellar populations that 
we are implementing in the semi-analytic models, but these 
data are not available to us. However we emphasize that the 
differences that arise from using the M05 models to convert 
from mass to light (or light to mass) are considerably larger 
than the ones origin ated from the determ ination of photo- 
metric redshifts (e.g. lMaraston et alj|2006h . The subsequent 
conversion from observed- to rest-frame magnitudes is more 
difficult to track, as the different theoretical templates usu- 
ally give differen t fitted ages depending on the properties o f 
each galaxy (e.g. lMaraston et al.|[200rj ; ICimatti et alJl200Sl) . 
However, this will produce differences between derived k+e 
corrections that are not systematic and therefore should not 
alter our conclusions. 

Finally, when comparing model results and data, one 
should consider that model magnitudes are "total", while 
observational measurements are usually based on "aper- 
ture" magnitudes. At redshift zero, a significant fraction 
of light might be missed for large o bjects that e xceed 
the available aperture di ameter (e.g. lLauer et all 120071 ; 
Ivon der Linden et alj I2007T ) . At the higher redshifts stud- 
ied here, despite galaxies being smaller than the maximum 
available apertures, there can still be an issue of missing 
light when small apertures are used to ensure high signal-to- 
noise. Moreover, the situation can be complicated by limited 
instrumental resolution that might blend together objects 
in crowded regions. The first problem is minimized in the 
data we use by applying poin t spread function (PS F) cor- 
rections to total magnitudes (|Cirasuolo et alj|20ld ). Nev- 
ertheless, both aspects can influence the evolution of the 
bright end of the luminosity function. 



4 RESULTS 

Fig. [2] compares the evolution of the 7f-band luminosity 
function from r edshift 3 to redshift 0.5 for the semi-analytic 
models with t helCirasuolo et "ail (l201df) d ata (shown as open 
black circles). |Pe Lucia fc BlaizotJ J2007I) models a re shown 
in red, MORGANA in green and lMenci et al.l (|2006h in blue. 
Original model versions are shown as dashed lines and the 
M05 versions as solid lines. 

The three galaxy formation models in the M05 versions 
show an enhanced if-band emission (between 0.25 and 0.5 
mags) from the brightest objects (Mk ^ —24) which for 
2^2 brings the models into agreement with data. The 
original versions of the models predict that only old popula- 
tions provide substantial 7f-band emission. For this reason, 
the bright end of the Tf-band luminosity function could only 
be built-up at lower redshifts, when old populations become 
dominant in massive galaxies. The TP-AGB phase gives a 
simple and straightforward way to solve the problem with 
the observed evolution of the AT-band. 

The agreement between the semi-analytic-plus-M05 
models for bright _R"-band objects and observations at 
high redshift is remarkable. In principle effects from the 
age/metallicity degeneracy could produce an artificial agree- 
ment between data and model. For example, a luminous 
Jf-band can originate from very metal-rich populations or 
from much older populations than those dominated by the 
TP-AGB emission. However, the wide redshift range that 



is spanned by these observations and the trend of the ob- 
served Tf-band LF with redshift allows us to exclude such 
effects acting at all redshifts. This is particularly the case 
at high redshift, where the time elapsed since the Big Bang 
is short enough such that this age degeneracy cannot enter 
the game. 

These results suggest that, if masses and ages were es- 
timated from observational data using M05, these would 
be in agreement with model predictions. The conver- 
sion to photometric properties was fully responsible for 
the dis agreement with obser vations that was pointed 
out by ICirasuolo et alj (I201C ) and previously found by 
other authors JPozzetti et al.l 120031 ; ICimatti et~all |2004| ; 
iKitzbichler fc Whitell2007l ). This result has important impli- 
cations for the observational determinations of stellar masses 
and ages from photometric data, in particular for galaxies at 
high redshift. Significant isT-band emission can be produced 
by young populations at high redshift through the TP-AGB 
stellar phase. Without considering it, large Jf-band emission 
can only originate at older ages, which results in a system- 
atic over-estimati on of stellar masses d erived from emission 
in this band (e.g. lMaraston et al.ll2006T ). 

Interesting diff erences among the mo dels emerge at low 
redshifts. The foe Lucia fc Blaizotl l|2007h model plus M05 
follows the bright tail in every redshift bin (z=0.5 to z=3.0). 
At redshift 1 the better match between data and the models 
implies a certain fraction of 1 Gyr population in these galax- 
ies, a prediction that could be tested by acquiring rest-frame 
near-IR spectra. 

For MORGANA and lMenci et al.l (|2006h . the inclusion of 
the M05 models worsens an existing discrepancy at lower 
redshift z <= 1.5, namely that the models over-predict the 
number density of massive galaxies. This discrepancy is em- 
phasized with M05 because existing ~ 1 Gyr populations 
have a higher flux. However, the excess is present for both 
versions of each model, hence it is primarily caused by mass 
growth rather than age. In both models, this is caused by a 
combination of enhanced merging times and inefficient AGN 

feedback at low redshift (see Section [2]). 

The dynamical friction merging times in M enci et al.l 

(2006) and MORGANA are shorter than what is expected from 
the numerical analysis of Bovlan-Kolchin et al.1 (|2008l ). as 
pointed out bv lPe Lucia et al.l ( 20101 '). Moreover, both mod- 
els have AGN feedback implementations that are inefficient 
in shutting down star formation in massive objects at later 
times. In the lMenci et al.l (|2006l ) model this is caused by the 
absence of "radio-mode" AGN feedback. In MORGANA some 
amount of star formation is required to destabilize the reser- 
voir of gas and trigger this feedback mode. This results in 
galaxy optical col ours that are too blue (as pointed out by 
iKimm et al.ll2009l . for MORGANA) and in an excess of massive 
galaxies as emphasised by both stellar population models in 
our study. 

In order to clarify the different results obtained with 
the three semi-analytic models, we show in Fig. [3]the model 
mass-weighted ages as a function of redshift for all galaxies 
brighter than M_r-=-24. At early times (z~2.5), the impact 
of the TP-AGB phase is larger, which stems from the mean 
ages around 1 Gyr of galaxies in the three models. As we 
move towards lower redshifts a bimodality emerges in the 
Munich model, with the bright-end of the AT-band lumi- 
nosity function being build up by a combination of young 
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M K (AB) Mk(AB) m k (ab) 

Figure 2. The evolution of the K-band luminosity function from z=3.0 to z=0.5. The original prediction s from three different sem i- 
analytic models are shown as dashed line s and their version with the M0 5 as solid lines. Red li nes refer to |Pe Lucia fe Blaizoti <2007ft . 
green to lMonaco et al". I J2007h and blue to Mcn ci et al. Data from Cirasuolo ct al. I ll2010h are shown as black opened circles and 

the solid line. 



and old populations (with the latter, maintained by the 
"radio-mode" AGN feedback, growing in importance as we 
move to low e r reds hifts). This bimodality is also present in 
iMenci et alj (|2006t ). but it is weaker and the oldest popula- 
tion is ~ 1 Gyr younger than its counter part in the Munich 
model. On the other hand, MORGANA shows considerably 
younger ages, centered at ^2.5 Gyr, with only a very week 
peak at ~ 4 Gyr. Despite the difference in the age distribu- 
tion at z—0.5, the ongoing star form ation produc e s the same 
mass excess in both MORGANA and lMenci et alj (2006) (re- 
sulting in a similar over-esti mation of the numb er density of 
bright iC-band objects). The|M enci et all (|2006t ) model has a 
smaller fraction of younger ages. This is due to the assumed 
modeling of "quasar-mode" feedback combined with the ab- 



sence of "radio-mode". At very high redshift (z>3) star for- 
mation is high in the progenitors of massive galaxies due to 
merger induced starbursts. At low redshifts only a fraction 
of these galaxies have their star formation quenched. On the 
other hand, MORGANA has continuous on-going star forma- 
tion but always at a moderate level, being self-regulated by 
the "radio-mode" feedback. 

It should be noted that the mass-weighted age is an av- 
erage over the individual populations that compose the theo- 
retical galaxies. This implies that individual ages can extend 
down to much lower values. This can be seen in Fig. [4] where 
we show the ages of the individual populations for galaxies 
with Mk < ~24 at z=0.5 in the three semi-analytic mod- 
els. The 1 Gyr old populations present in the three models 



© 2009 RAS, MNRAS 000, [T]{9] 



The rest-frame K-Band Luminosity Function 7 



0.08 
0.06 
o 0.04 



Moss-Weighted Ages 
M„<-24 



0.25<z<0.75 




Age (Gyr) 



10 



Figure 3. The mass-weighted age distribution of galaxies with 
M-K < —24 in the three semi-analytic models used in our study. In 
the th ree panels, the solid red line represents |Pe Lucia fc B laizot 
(2007), the solid green line shows MORGANA p r edicti ons and the 
solid blue line gives the ages in iMenci et al,l | |2006| ) . From top 
to bottom, the panels show respectively galaxies ages at z=0.5, 
z=1.0 and z=2.5. 



explain the impact of the TP-AGB phase even at this red- 
shift. It can also be seen that MORGANA exhibits the larger 
fraction of young stars, since the "radio-m o de" f eedback is 
regulated by star formation. IMenci et alj (I2006I ) ages are 
considerably older, due to the impact of the "quasar-mode" 
feedback at early times. However, as it starts being inef- 
fective at lower redshifts a considerable fraction of younger 
populations emerge. 

Despite the significant improvement obtained in match- 
ing the evolution of the bright end of the i^-band luminos- 
ity function, the faint end remains problematic. The lumi- 
nosity (as well as the stellar mass) function for faint ob- 
jects (—22 ^ Mk ^ —24) is known to be much higher 
than measured jFontana et a l. 2006; Weinmann et al. 2006; 



iHenricmes et al.M200Sl ; IFontanot et all l2009l b) and the 
elusion of the M05 models worsens the case. This excess 
can be removed in different way s at redshift zero, by us- 
ing a more up-to-date cosmology (jSomerville et al1l2008l ) or 
combining the d isruption of stellar material from satellites 
during mergers (M onaco et al.ll2006l ; iHenriques et~ai] 2008; 
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Figure 4. The ages of the individual populations present in semi- 
analytic galaxies wi t h My < —24. The solid red line represents 
|Pe Lucia fc Blaizoll ((2007), the solid gr een line shows MOR GANA 
predictions and the solid blue line gives IMenci et al ] l|2006h ages 
at z=0.5. 



ISomerville et aL 2008f) with more efficient superno va feed- 
back (|Henriques fc Thomas! l20ld ; iGuo et al.ll2010l ). Never- 
theless, the comparison pres ented in our work for high red- 
shift (as already shown by IFontanot et all 12009b ). shows 
that for the early phases of galaxy evolution this might be a 
problem, even considering problems of incompleteness with 
the high redshift data. 



5 SUMMARY 

The main objective of this work is to re-address the fun- 
damental question of matching the observed rest-frame K- 
band luminosity function of galaxies over th e Hubble time, 
using s emi-analytic m o dels. I n the literature l|Pozzetti et alj 
|2003l; ICimatti et al l |2004| ; iKitzbichler fc White! 120071 ; 
ICirasuolo et alJ l2010l ) , it has been pointed out that semi- 
analytic models underestimate the rest-frame A"-band 
galaxy luminosity of the brightest objects at high redshift 
(~ 2 — 3), and the failure has been attributed to an insuf- 
ficient mass build-up at early epochs. 

However, the galaxy luminosity function does not only 
depend on the mass build-up, but also on the light emit- 
ted per unit mass. Hence, in order to pin down the origin 
of the mismatch, we improve upon the rest-frame if-band 
emission from the model galaxies. We use the M05 stellar 
population models, which include emission from the cool and 
luminous TP-AGB phase. The contribution of this phase of 
stellar evolution in the M05 models is important at inter- 
mediate ages (between 0.2 and 2 Gyrs), which are expected 
to be predominant at 2<z<3. The relevance of this ingre- 
dient has b een recently shown f or th e semi-analytic model 
GALICS in iTonini et~aTJ l|2009l . |2010| ). where the observed 
near-IR colours of redshift 2 galaxies could only be matched 
by the model inclusive of t he TP-AGB emission. Similarly, 
IFontanot fc Monacol (|2010f ) showed that the inclusion of this 
stellar phase in MORGANA increases significantly the num- 
ber density of EROs at high redshift. 

We consider several semi-analytic models - namely 
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\Be Lucia fc BlaizotJ (120071 ') , MORGANA (Monaco et al. 2007) 
and Menci et al. ( 20061 ) and we implement the M05 stel- 
lar population models, keeping all other ingredients and as- 
sumptions unchanged. 

We find that the semi-analytic models with the M05 
models exhibit a brighter K-bsmd LF by as much as 0.5 
mags at the highest redshift bins. This is precisely the offset 
that was plaguing the com parison with the UKID SS data 
for the brightest objects in lCirasuolo et all (|2010h . Models 
and data at high redshift and for Mk < —24 now match 
very well. 

This result is strongly suggestive that the models at 
redshift 2 — 3 do not underestimate mass, rather they did 
require a proper conversion between mass and light. More- 
over, we show that a strong evolution in mass predicted by 
hierarchical models is compatible with no evolution on the 
bright-end of the Jf-band LF from z=3 to the local universe. 
This means that, at high redshifts and contrary to what is 
commonly accepted, A"-band emission is not necessarily a 
good tracer of galaxy mass. 

At lower redshift, the details of the implementation of 
AGN feedback and merging timescales produce differences 
between the various semi-analytic models that are not al- 
tered by t he inclusion of the M 05 models. In particular, MOR- 
GANA and Me nci et al.l (|2006i ) exhibit a if -band luminosity 
function bright tail that is higher than the data, which is due 
to an excessive mass build-up connected to the lack of an 
efficient quenching of low-z cooling flows via "radio- mode" 
feedback. 

Similarly, the faint end of the galaxy luminos- 
ity function remains substantially overestimated by 
the models at a ll redshifts. This is a well docu- 
ment e d problem llFontana et aT 20061; Weinmann et al.l 



2006; Henriques ct al. I 120081: iHenriques fc Thomas! l201oT 
Fontanot et al.l 120091 b: iGuo et al.l l2010h that we plan to 
study in future work. 

In recent years, our understanding of the various phases 
of stellar evolution has improved. Moreover, we now have 
high-quality observational data covering a wide spectral 
range (including the rest-frame near infra-red). Therefore, 
we should now be able to constrain galaxy formation mod- 
els with better accuracy and disentangle between different 
theoretical approaches. 
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